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Abstract. 

By demanding that a bounce is nonsingular and that perturbations are well-behaved 
at all times, we narrow the scope of possible models with one degree of freedom that 
can describe a bounce in the absence of spatial curvature. We compute the general 
properties of the transfer matrix of perturbations through the bounce, and show that 
spectral distortions of the Bardeen potential <i> are generically produced only for the 
small wavelengths, although the spectrum of long wavelength curvature perturbations 
produced in a contracting phase gets propagated unaffected through such a bounce. 
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1. Introduction 

It has become generally admitted, especially with the recent WMAP data [1], that 
the Universe must have undergone a phase of inflation |2], i.e. a very short period of 
time during which the ongoing expansion was exponentially accelerated. This phase 
not only solves the usual cosmological flatness, homogeneity, monopole excess and 
horizon problems, but it also produces, as a bonus, an almost scale-invariant (usually, 
but not always, slightly red) spectrum of primordial scalar perturbations. These small 
inhomogeneities, of one part in about 10^, have the right spectrum and can be given, 
with some amount of fine-tuning, the right amplitude to seed the large scale structure 
formation. In fact, it can be argued [3] that the inflationary expansion and the 
ensuing superadiabatic amplification of the zero-point energy of the quantum fields 
is the only plausible mechanism to transfer microscopic quantum fluctuations up to 
the cosmologically relevant length scales without annihilating the amplitudes of those 
fluctuations in an ever expanding universe. 

Bouncing models [U [5] have been proposed as alternatives to this scenario, mostly 
in the framework of string theory [6],[7j (see, however, Ref. [8j). A bounce, i.e. a period 
of contraction followed by expansion, could explain the flatness if the expansion phase 
lasted much less than the contracting era; it could explain homogeneity by making the 
past light cone very large during the contracting era so thermalization could take place; 
and, as it turns out, it can very easily give rise to the same mechanism of superadiabatic 
amplification as inflation [9j. 

The main distinction of bouncing models compared to inflation lies in which 
term dominates the spacetime curvature R = —12H^ — 6H. Whereas in inflation 
H <^ and the physical wavelengths grow much faster than the curvature radius 
~ p ' , close to a bounce if ~ and the curvature radius grows to 
~ H^^^^ ~ [—{p + p)]~^^'^ as the contraction rate grinds to a halt, then falls back 
down rapidly as the expansion phase begins. 

This means that the modes of interest are pushed inside the curvature radius during 
the bounce, and then out again as the universe expands. This "in-out" transition is what 
makes superadiabatic amplification possible, both in inflation as well as in bouncing 
models. The question is whether sufficiently natural models can be found which give 
rise to near-scale invariant spectra of cosmological perturbations [9]. 

As opposed to inflation, in which the phenomenological consequences of the 
simplest single-field, slow-roll models are extremely similar (slightly red spectra), 
in bouncing models the ensuing spectrum of cosmological perturbations can vary 
dramatically, depending on the model. Moreover, making the universe bounce is far from 
straightforward since general relativity forbids this behavior as long as the Null Energy 
Condition (NEC) holds. As a result, bouncing models can become rather intricate. The 
simplest bouncing models developed so far have relied on a combination of fluids |10j . 
the presence of spatial curvature [12j or ghost fields [13j. Some models predict mode 
mixing with or without spectral modifications through the bounce itself, and it has been 
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suggested that these features essentially originate from either the mixture of two fluids, 
i.e. from the entropy perturbations (see ^0] and, in particular, [11] where the precise 
treatment of entropy perturbations is done), or from spatial curvature [12]. Hence the 
need for a single field flat space bouncing model. 

Therefore, we propose here a minimalistic model for a bounce with a single matter 
component (a generalized scalar field, or K-essence) and zero spatial curvature. The 
condition that spatial curvature is small shortly after the bounce is a natural one if 
the contracting phase lasted much longer than the ongoing expansion phase. We also 
demand that the energy density is positive at all times, that a non-singular bounce takes 
place, and that the sound speed of perturbations is well-behaved at all times. 

By expressing the Lagrangian of the generalized scalar field as a Taylor series around 
the field and its momenta, we can easily implement these constraints and proceed 
to construct a very general class of sensible bouncing models with a single fluid and 
no spatial curvature. The only shortcoming of our class of models is that, because 
H > near the bounce, they all lie in the "phantom" sector, w = p/ p < —1, so the 
connection with an expanding radiation era would necessitate the introduction of matter 
fields and a decay mechanism similar to preheating [H]. This is precisely the scenario 
recently proposed in Ref. [15], where an explicit scenario is realized using the ghost 
condensate model [16]. For additional context on the use of non-canonical scalar fields 
in cosmology, see, e.g., Ref. [171 [THl HQ] . Note also that if one assumes a contracting phase 
dominated by normal matter (preferably pressureless matter, in order to get an almost 
scale invariant spectrum of perturbations [20]), then because the phantom divide cannot 
so easily be crossed [21], there must also exist a transition between this contraction and 
our K-bounce, equivalent to preheating but in the other way, that one could henceforth 
call precooling. 

The advantage of our models lies in the simplicity of their perturbative sector. We 
show explicitly that cosmological perturbations can be propagated in a non-singular way 
through the bounce. We also show that the perturbations are well-behaved through the 
numerous instantaneous de Sitter phases (moments of time at which H = 0) that take 
place in our model. 

We have computed the transfer function for perturbations, and we show that an 
initial spectrum of cosmological perturbations can get distorted by the bounce. As 
this distortion depends on the duration of the bounce, our conclusion is that bouncing 
models generate power spectra with a wide variety of scale dependences. However, the 
scale dependence of the transfer matrix is important only for short wavelengths, so that 
the cosmologically relevant (large) scales are transferred through the bounce unaffected. 

This paper is organized as follows. First, we describe the K-essence model 
generating the K-bounce, provide the relevant equations of motion for the background 
and derive the conditions under which a bounce is possible (Sec. II). We then specify, 
in Sec. HI, through a Taylor expansion around the bounce, the form of the pressure 
function we use afterwards, and provide the constraints for a non-singular bounce to 
take place. Sec. IV discusses a number of specific background models and attempts at 
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classifying them by the bounce duration. We then move on, starting in Sec. V, to the 
study of the perturbations. We first reduce the overall system to a single equation for 
the only degree of freedom avalaible, which we chose to be the Bardeen gravitational 
potential We discuss analytically several potentially problematic cases (the bounce 
itself, the instantaneous de Sitter phase, and the quasi de Sitter bounce), and we show 
that $ is well behaved at all times. Having shown the propagation of linear perturbations 
across the bounce to be regular at all times, we then compute numerically this time 
evolution, setting initial conditions at an arbitrary time at which we impose $fc = 1 
and = for the Fourier modes. We end up with some considerations about model 
building in a concluding section. 

2. Generalized scalar-field models 

We will assume that the matter sector is represented by a scalar field Lagrangian of the 
form 

£ = v^p(X,0), (1) 

where 

X = ^ g>^^d,(t^ , (2) 

and use the timelike signature, diag(+, — , — , — ) for the metric g. From this Lagrangian, 
one gets a stress-energy tensor reading 

T^"' = {p + p) u^'u" - pg"", (3) 

with energy density p = 2Xpx — P and = 0^^/ V2X. These relations give back 
the usual one for the canonical scalar field theory provided one then takes the simplest 
Lagrangian function = ^ — ^(0)- 

Since Ostrogradski's theorem [221 |23l [21] precludes local higher derivative terms 
from appearing in the action principle, Eq. ([1]) is the most generic scalar field Lagrangian 
which may be stable. Notice that our Lagrangian does not need to be separable in 
terms of functions of the kinetic term X and the field (p, as is sometimes assumed for 
K- inflation [2S] or K-essence [27] . 

The important aspect of the quantum instability of this theory would also need 
to be addressed, since evidently any Hamiltonian which is unbounded from below 
would be instantly destroyed by quantum tunelling of positive-energy particles into 
the negative-energy particles [21]. For theories with non-canonical kinetic terms the 
quantum stability is a nontrivial issue, in particular for the case of "phantom" models 
- see, for instance, Ref. [21] . 

Introducing the fiat Friedman-Lemaitre-Robertson- Walker metric with scale factor 
a{t), namely 

ds^ = dt^ - a'ifydx^, 
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and the Hubble parameter H = a/a {a. dot standing for a derivative w.r.t. the time 
coordinate t), the Einstein field equations are then given by 

3H^ = SttG p = 8ttG {2Xp^x - p) , (4) 

-3H^ - 2H = 8nG p . (5) 

As for the matter field, the Euler- Lagrange equation stemming from Lagrangian ([1]) is 
nothing but the conservation of the stress-energy tensor ([3]), i.e. 

V^T^^ = =^ p,xn<f) + d^d^{p,x)-p,^ = 0, (6) 

which, under the assumption of homogeneity of the scalar field, 4>(t), is reduced to 
the simpler form, 

^ + 3if0+M=O, (7) 

which reduces to the Klein-Gordon equation for the canonical theory Cc- Here the sound 
speed is given by 

4 = ^ = — ^-^ , (8) 

^ p,x 2Xp^xx+P,x ^ ' 

and it should be clear from the unapproximated equation of motion that it is the function 
responsible for the speed with which inhomogeneous scalar field fiuctuations propagate 
through spacetime. In particular, a negative Cx would give rise to exponentially growing 
small-scale fiuctuations, meaning that the theory is classically unstable. 

Since our final aim concerns the predictability and spectrum of cosmological 
perturbations before and after the bounce in a one-fiuid model, our first requirement is 
that the sound speed never becomes negative. We also demand that it remains finite, 
since a diverging sound speed would cause a singularity in the transfer matrix [121 US] 
that relates the cosmological perturbations before and after the bounce - destroying, 
once again, the predictability of the theory. Therefore, our first physical constraint is 

< 4 < oo . (9) 

Notice that even though we demand that the sound speed squared is always positive 
and finite, we should still work under the assumption that our models are just 
phenomenological realizations of some unknown fundamental theory, so that the second- 
quantized perturbations of the gravitational degrees of freedom are not being taken 
into account properly here. Otherwise, since both p^x and 2Xp^xx + P,x are negative 
through the bounce phase in our models, the theory can become unstable, decaying 
instantaneously through graviton production [2^ [28]. 

Our second requirement is that the energy is non-negative. In particular, if the 
bounce happens at t = we must have that 

p(to) = 0, (10) 

and, as a result of the Einstein equations written in the form — 3i/^ — 2H = SnGp, we 
conclude that we must impose 

p{to)=po<0. (11) 
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This means that, ast ^ t^, the equation of state parameter w = p/ p becomes infinitely 
negative. But w < —1 is the domain of the so-called "phantom" (or "ghost") models, 
and it has been shown that crossing the "A barrier" w = — 1 is impossible in simple 
single-field models [2T| [29] . Therefore, our bouncing model is limited to w < — 1, 
which in practice means that in the asymptotic past (future) the Universe approaches 
a contracting (expanding) de Sitter stage. These limiting stages must somehow be 
connected with non-phantom dominated epochs through precooling and preheating 
phases. 



3. Taylor expansion of the Lagrangian 

Our model relies on a series expansion of the Lagrangian in terms of the field and its 
momentum. Without loss of generality we set the value of the field at the bounce to be 
(p(to) = 0, and its time derivative X(to) = Xq ^ 0, so we can write 

p{X, (f))=po+ Px {X -Xo)+p^(j) + PX4, <P{X - Xo) (12) 

+ \pxx {X - XqY + ^p^^ 0' + • • • 

In order to obtain a well-behaved bounce, it is helpful to assume that the behavior of 
the field near the bounce is analytic in time 

0(t) ^ 01 (t - to) + 02(t - t^f + 03(t - tof + ■ ■ ■ , (13) 

which means that Xq = |0^. We stress that the Taylor expansion in Eq. (fT3|) is not 
used in any way to constrain the dynamics - we only use it as a means to adjust the 
parameters of the Lagrangian in light of the constraints. 
The constraint that p(to) = translates into 

2pxXo -po = =^ Po= Px<Pl ■ (14) 

A second constraint comes from the stress tensor conservation, i.e. p = —3H{p + p), 
imposing that p — > at the bounce. This means that only the term in p which is 
quadratic in time survives. In terms of our parameters, this condition is expressed as 

p^ = 202PX + 0?(PX0 + 202Pxx) • (15) 

Finally, the Friedman equation SvrGp = —3H^ — 2H at t = to leads to the third 
constraint, namely that px < and po = px<Pi < 0, with px given by 



Px 
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Notice that 0i and 03 must be chosen such that the square root is real, and such that 
Px is negative. Notice also that the + and — branches are identified by simultaneously 
changing the signs of 0i and 03. 

To summarize: our set of constraints determines some relationships between the 
Lagrangian parameters Pq, p^ and px in the context of the class of models in which 
the behavior of the scalar field near the bounce can be represented as Taylor series. 
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Presumably, other models for which the scalar field around the bounce cannot be 
represented by such a series will lead to similar constraints between the parameters 
involved in these cases. As we are interested in some minimalistic bouncing model and 
the general conclusions that can be drawn thereof, this will suffice for us. 

Notice that the scalar field parameters (pi are essentially free, and that the 
Lagrangian parameters p^^f), px<j> and pxx are also essentially free - we only need to 
make sure that the square root in Eq. f|T6|) remains real. Higher-order parameters 
(yPxxx etc.) would come into these constraints, but they would also remain basically 
free. This means we can tune the parameters of the Lagrangian in order to make the 
models stable - which is very important for phantom models. It also means that we can 
set the parameters such that the bounce is short or long, fast or slow, at will. 

4. Concrete models: background 

The class of models one can construct with the procedure above has a very rich 
phenomenology. In all of them the conditions we impose on the parameters are such that 
any bounce, defined as the point in time at which H = 0, is necessarily non-singular, 
having if > at this point. Therefore, if we set our initial conditions to a Universe 
that is contracting, it necessarily will end up bouncing provided the constraints on the 
underlying parameters are indeed satisfied. 

There is only one kind of fixed point in our theory, namely H ^ 0. As a result, 
and whatever the initial conditions, once we have passed through the bounce, our 
models necessarily asymptote to a de Sitter Universe; this fixed point is an attractor 
provided H > 0, and a repulsor otherwise. In practice, some intermediate quasi-de Sitter 
phases (contracting as well as expanding) can happen as the model contracts and then 
expands, which is rather interesting from the point of view of the background model, 
but represents a formidable complicating task if one is interested in the perturbations. 

We have chosen to concentrate on three concrete models, one in which the bounce 
is relatively fast and short, one in which it is a slow and long phase, and another in 
which we tuned the parameters so that the bounce is also a quasi-de Sitter phase (i.e., 
both H = and if = at the bounce.) All models approach a contracting (expanding) 
de Sitter phase in the past (future), which is natural since going backwards in time 
transforms the repulsor with H < into an attractor. The contracting phase is in 
fact an unstable point which all trajectories exit from, whereas the expanding de Sitter 
phase is an attractor point where all our models must finish. Therefore, in order to make 
the transition to a radiation-dominated Universe we must introduce new ingredients, or 
make the scalar field decay into some other fields. As this reheating process usually 
preserves the basic properties of the cosmological perturbations (at least in the single 
field case at hand), we will not treat it here - see, for instance, Ref. [15]. Similarly, if 
we want to originate with a stable-matter dominated phase, we will need a transition 
{precooling) to lead into the bounce phase. For the same reasons as the preheating, we 
shall not consider the details of such a transition, and will just assume, as usual, that 
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the long-wavelength spectrum of perturbations is transmitted unchanged through these 
precooling and preheating phases. Technically, this translates into saying that we do 
not impose physically motivated initial conditions here, assuming that they have been 
generated in the phase preceeding the precooling, and therefore out of the scope of this 
paper. 

It is useful to write down a few identities for the background that hold in general. 
The equation of state can be written, with the help of Friedmann equations, as 

w = - = -l--—, 17 
p 3 

while its time derivative can be conveniently expressed as 

w H H 

-'^-Tr^ (18) 



1 + w H H 
where we have used the continuity equation, p = —3Hp{l + w). 

5. Perturbations 

We now perturb the scalar field as <^ ^ (j){t) + 6(j){x,t), and for the metric we fix the 
gauge to the conformal-newtonian (longitudinal) one as [25j 

ds^ = [1 + 2$ {x, t)] dt^ - [1 - 2$ {x, t)] a^{t)dx^ . (19) 

By using the constraint equations in the case of a single generalized scalar field we can 
express the Mukhanov- Sasaki variable [26] as 

v = zC, (20) 

where 

c = $ + = $ - 5 ($ + //$) , (21) 

and in our case (w < —1) we have 



2 c-^ 

The Mukhanov variable v obeys the equation 

v" + (^cj,e -^^v = 0, (23) 

where a prime denotes a derivative with respect to conformal time 1] = J dt/a, 
i.e. d/dt] = ad/dt. In lieu of Eq. f l23|) we can view z" /z as an effective potential that is 
scattered by the incoming wave v. 

However, it can be immediately seen from Eq. (!22|) that the transformation to the 
variable v is ill-defined in two particularly important situations: first, if the equation 
of state goes to infinity, as happens in our bounce, and second, if w — > — 1, as happens 
if the Universe reaches a de Sitter phase. In fact, the effective potential z"/z becomes 
singular in these situations. Obviously, in these cases the Mukhanov variable cannot be 
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usefully employed and we must search for other, more suitable ways to represent the 
perturbations. It is interesting to realize that the situation is similar to what happens in 
the curvature dominated bounce examined in Ref. [12]: the Mukhanov variable becomes 
useless in both these bounce cases. We must therefore resort to the original Einstein 
equations for the metric perturbations directly. 

In the single scalar field case it is possible to write a second-order differential 
equation for $, which reads 

*-.(^-f)i-H(4^ + 2^f-f^)*^0. (24) 

It is clear that this equation is completely well-behaved through a bounce (i/ — 0), as 
long as H remains finite. 

As is well known, this equation also describes well the perturbations in a nearly 
de Sitter (inflationary) spacetime. This is evident if we write Eq. (|2^ in terms of the 
slow-roll parameters e = —H/H^ and 5 = —H/ (2HH) 

$ + (1 + 25) + [c\-^ - 2e + 25^ iJ^^ = . (25) 

It then becomes obvious that by taking e and and neglecting the 

exponentially decaying gradient term we obtain that in the slow-roll regime 

<l>s.R. ~ Ae-(^-^+2e)m ^ B^-2(5-e)m ^ (26) 

to first order in the slow-roll parameters. 

However, Eq. (1211) may not be appropriate in an instantaneous de Sitter point, 
i.e., an instant of time when w = —1 {H = 0.). We have found that it is particularly 
enlightening to write the following set of first-order equations 

=4^*. (27) 
t + i/4= -^(C-*). (28) 



These relations, together with Eq. (1211) . show that the perturbations are propagated 
through the many important phases described below in a regular way. 

5.1. Exact solution near the bounce 

First, consider a bounce at t = that occurs within the class of models given by Eqs. 
([I2D-([T3]13. We can then write 

H^H,t + iH^t" + ^ifst' + ■ ■ ■ , (29) 
2 D 

I Here and in the following subsection, the choice t = for the point under consideration is of course 
a mere convention aimed at simplifying the subsequent equations; in the numerical approach, we will 
set the initial contracting solution at < = 0, so the bounce takes place at a different location. 
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where we take iiTi > in accordance with our class of models for which 1 + w < 0. 
We suppose the approximation above to be valid for small times such that \t\ <ti 
min {2Hi/\H2\,3\H3\/\H2\). 

Substituting the approximation above into Eqs. (127|) - (!28|) or, equivalently, into Eq. 
( IMl) . and keeping only the dominant terms we obtain the following equation for $: 

^+ + 2i/i + H^t] $ ^ , (30) 

where is the (near-constant) scale factor at the bounce, so we implicitly assume that 
\t\ <^ H-^^\ We will also assume that the sound speed Cx is approximately constant 
during the bounce, which is the case in all models we have considered. It is interesting 
to notice that if H2 = 0, then there is a limiting case H3 = Hf {H = HH) for which 
the exact solutions near the bounce are pure oscillatory modes. 

By writing the solution to Eq. fl30l) as a truncated Taylor series in time, it is easy 
to find two linearly independent approximate solutions, 

*i « 1 - (1 + 111) H,e- - ^ (3 + T?) , (31) 

where 7^ = c^k"^ /{alHi). We can in fact find exact solutions to Eq. (!3Ql) . and the 
two linearly independent modes turn out the be essentially a Hermite polynomial and a 
confluent hypergeometric function $. Both functions are analytic at t = and reduce, 
to lowest order in t, to the approximate solutions fl^ . 

In terms of the curvature fluctuation the approximate solutions are: 

Ci - 1 + l^lH.f - flit' , (32) 

Therefore, the curvature perturbation is completely regular across the bounce. Notice 
that neither the growing nor the decaying modes of the curvature fluctuations near the 
bounce depend on the cubic term in Eq. fl2U]) . even though the newtonian potential 
$ does. From Eqs. (!2T!) and (13T]) - (!32|) we can also see that by keeping only the dominant 
mode we make the curvature fluctuation ( equal to $ at the bounce. 

Notice that what was the growing mode in the contracting era becomes the decaying 
mode in the expanding era, and vice-versa. This behavior is completely generic for the 
linear cosmological perturbations, and has been shown to work in much more complex 
bouncing models [30] . 

5.2. Exact solution near an instantaneous de Sitter phase 

Now we analyse the solution near a de Sitter point - i.e., and instant of time when 
w = —1. Since we assume that the sign oi 1 + w = —2H/3H'^ does not change, we take 



1 + 



Hi 
Hi 



El 

Hi 



H2 
Hi. 
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the following approximation for the Hubble parameter near a de Sitter point which we 
place at t = 0: 

H^H, + ^H^t^ + ^H,t' + ■■■ , (33) 

where we suppressed the linear term because we want H = at t = 0, and the 
absence of a quadratic term is implied by 1 + w < 0. The approximation is valid 
fo^ \M ^ \^Hq/ Hs\^^^; we have found that it is always valid near de Sitter points in our 
class of models. 

Neglecting the subdominant terms we obtain the following equation for $: 

^+(^A-j^^+(^B- $ = , (34) 

with 

A = Ho - ——, and Bk = 



3H3' 3H3 ■ 

Notice that the cubic term H3 does not show in the equation for $ at leading order 
order - it will, however, reappear when we compute (. 
Defining 



fifc = ^l-4^ ' ^^^^ 

rescaling the time to z = AQkt, and making the variable change $ = 
exp [—(1 + VLk)z/ (2f2jt)]?/(z), we can reduce Eq. (134|) to the equation for the confluent 
hypergeometric function, 

zy" + (7 - z)y' - ay = , (36) 

where 7 = 4 and a = 2 + {2Hq — A)/{AVtk). The two linearly independent solutions to 
Eq. (136!) are given by 

yi = iFi(a,7,2;) , (37) 

00 

y2 = C iFi {a, 7, z)\io.z + z~^ ^ t;„z" , 

n=0 

where the coefficients C and u„ can be found in standard textbooks on special 
functions [31]. Notice that the first solution is well-behaved everywhere, but the second 
solution is non-analytic at the origin due to the presence of the log term. This happens 
because the confiuent hypergeometric function iFi(a,7, z) with integer 7 has a branch 
cut in the Riemann plane. As a result, we could not have found the second solution 
by writing a naive Taylor series around t = 0, as was done in the previous section. 
Nevertheless, specifying the value of the second function and its derivative anywhere 
fixes its value everywhere in the Riemann plane, so the solution can be propagated from 
negative to positive values of t. Of course, Eq. flM|) is only approximate, so the exact 
solution to the exact equation may be much better behaved, but this subtlety rendered 
the actual numerical evolution of the perturbative equations tremendously complicated 
at de Sitter points. 



K-Bounce 



12 



In terms of the newtonian potential we find tlie following approximate solutions for 
small z, namely 

$1 ^ (4 - z) + 0{z^) , (38) 
$2 ^ 1 - -s; + C2Z^ - csz^ + (c4 + dilnz) z"^ -\ 

The first is a decaying (growing) mode in the contracting (expanding) phase, while the 
second is a constant mode. 

Upon substitution of the modes ( l38i) into the definition of ( we find 

Ci ^ ai - bit^ + 0{t^) , (39) 
(2^t-'[a2-b2t' + 0{t')] . 

The first solution is nothing but the constant curvature mode that passes essentially 
unaltered through the instantaneous de Sitter phase, while the second solution has a 
pole ~ at the de Sitter point but has no constant piece. This pole corresponds to 
no real physical singularly: it just points out the inadequacy of the definition of the 
curvature fluctuation in this situation since the solution for the newtonian potential $ 
is completely well-behaved and can be propagated through any de Sitter point. 

5.3. Exact solution near a quasi-de Sitter bounce 

An interesting limiting case happens when the bounce occurs in such a way that, as 
if 0, — ^> as well. Since if > in our types of models, we conclude that near this 
quasi-de Sitter bounce we also have H ^ 0. Therefore, near the quasi-de Sitter bounce 
the Hubble expansion parameter can be expanded as 

H ^ ]:H,t^ , (40) 

which, substituted into Eq. fl24l) . leads to 

^-h + (3l<i>^0, (41) 

with Pi = Cxk'^/al. There are trivial solutions to this equation in terms of spherical 
Bessel functions j_|_3, 

, /COST . \ ^ f sinr\ 
^ ^ At y h smrj + 5r ( COST j , (42) 

where r = Pkt- The approximate solutions for the curvature perturbation ( around the 
bounce are then given by 

C . ^ (l + Ir') + B (-1.=) . (43) 

So, again we see the presence of a constant mode, and of another mode which decays 
rapidly in the contracting phase but grows rapidly in the expanding phase. 
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5.4- Numerical evolution in some concrete models 

We now study a few concrete models. It is important to keep in mind that the condition 
that 1 + w < (or, equivalently, > 0) at the bounce means that the equation of state 
must be phantomhke at all times, as crossing the phantom barrier is prohibited in single- 
field models [21]. Since in the phantom case the only stable fixed point is w — > — 1, 
it is only natural that the asymptotic solutions are quasi-de Sitter (contracting and 
expanding), as we indeed find. 

Nevertheless, our freedom to set the parameters of the Lagrangian (po, P(/>, Px, etc.) 
means that we can vary the duration of the bounce. Since if = at the bounce, the 
time scale which sets how fast the bounce occurs is naturally given by ATb ~ a/ Hb- 
Therefore, by tweaking Hb we can construct models in which the bounce is very fast or 
very slow. 

As we are ultimately interested in the cosmological perturbations and their spectra 
in these models, it is useful to consider what should happen to the perturbations if 
the bounce is very fast or if it is very slow. Eq. ( |23i) tells us that we can regard the 
problem of the propagation of cosmological perturbations as that of the scattering of a 
wave function f by a potential = —z"/z. By changing the duration of the bounce, 
we are in effect changing the potential and changing the interval of time in which 
the wave function interacts with the potential. Hence, intuitively we should expect 
that for very fast bounces the short wavelengths will barely reach reach the oscillatory 
regime, while for slow bounces the oscillating stage will be fully realized by the short 
wavelengths. Hence, we should expect that, for those wavelengths that can reach the 
oscillatory regime, the change in their amplitudes is going to be more drastic in slow 
bounce models than in fast bounce models - see, later. Fig. 5. 

We have constructed three models which are broadly representative of the 
phenomenology of K- matter bounces: a fast bounce (FB), a medium bounce (MB) 
and a slow bounce (SB). 

First, consider the fast bounce (FB) of Fig. 1. As shown in the upper panel, the 
universe starts in a quasi-de Sitter contracting phase, with w = —1 and H ~ —1.1. 
It contracts with that initial rate up until t ~ 4.8, then it bounces at t ~ 5 as the 
expansion rate grows very rapidly. It then reaches another quasi-de Sitter phase, albeit 
an expanding one. In our arbitrary time units, the bounce lasts about At ~ 1. Notice 
that at t ~ 6 the density and the expansion rate become flat for an instant of time, 
meaning that at that point the equation of state reached the value w = —1 - that is, 
the universe went through an instantaneous de Sitter point. Notice also that nothing 
special happens to the sound speed c\ - indeed, in all our models the sound speed is 
well-behaved and is not crucial to any of our discussions. 

The cosmological perturbations in the FB model are shown in the lower panel of 
Fig. 1, for two different wavelengths {ki = 10~^ and k2 = 10), along with the scale 
factor. As discussed above, we do not have a natural criterium to impose on the initial 
conditions of the Bardeen potential. Thus we chose, for all numerically evolved models 
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Figure 1. Fast bounce model (FB.) In this model the parameters are: po — —3.815, 
P0 = -0.891, px = -4.710, = 1.0, P4,x = -2.4, pxx = 0.5 (with /i = 0.9, 
/2 = 0.1035 and /a = —1.) Upper panel: background quantities H, p, 0, and 
p + p. Lower panel: scale factor a and perturbations in the metric ($) for two different 
wavelengths, fci = 10^^ and k2 = 10. Notice the kinks in the absolute value of ^fc^, 
which indicate a regime of oscillations. The kinks in , on the other hand indicate 
that it starts positive, then, during the contracting era at around t ~ 2.5 it becomes 
negative, and then it becomes positive again in the expanding era at t ^ 11. The kinks 
in $fcj are a manifestation of the changing roles of the decaying and growing modes 
before and after the bounce. 



below, to set •I'ini = 1 and ^ini = 0: then, getting anything else but a constant $ for 
asymptotically long times after the bounce would be evidence of mode mixing. 

It can be seen that the newtonian potential is well-behaved at all times (as shown 
in the analytical solutions of the previous sections). For wavelengths longer than that 
of the mode ki the solutions for the perturbations are all identical, meaning that the 
bounce does not affect them differently. With our initial conditions all perturbations go 
through a sign change at around t ~ 2.5, which is just a manifestation of the relative 
growth of the dominant mode compared to an initially mixed-mode state. 

Notice that the small- wavelength mode k2 detaches from the behavior of the mode 
ki at around t ~ 5, which indicates that the small- wavelength mode almost reaches 
the oscillatory regime. Indeed, for wavelengths smaller than that of the mode the 
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t 

Figure 2. Medium bounce model (MB). In tliis model the parameters are: po = 
-1.776, p0 = -1.915, px = -1.776, p^^ = 2.3, p^x = 0, pxx = 0.5 (with /i -1, 
/2 = 0.75 and /s = 1.) Upper panel: background quantities -ff, p, 0, and p + p. 
Lower panel: scale factor a and perturbations in the metric (<&) for two different 
wavelengths, fci = 10^^ and k2 ~ 10. 

perturbations go through a period of oscillations which becomes longer as we consider 
smaller wavelengths. This means that these modes are small enough to be insensitive 
to the curvature radius created by the bounce. Equivalently, we can say that the modes 
V experience a very small effective potential Vy, so they simply oscillate. 

In Fig. 2 we show a bounce model (MB) in which the bounce itself happens over 
a longer period of time, At ~ 4. We have also set the parameters so that the instant 
of the bounce coincides with an instant when p + p — > 0. Hence, in this model the 
bounce {H = 0) is also a quasi-de Sitter point (iJ = 0); in other words, the background 
behaves, close to the bounce, like Minkowski spacetime. It is interesting, although not 
entirely unexpected, that even in this critical model the perturbations are entirely well 
behaved at all times. 

In Fig. 3 we show the SB model. Here the bounce is accompanied by many quasi-de 
Sitter instantaneous points {p + p = 0.) The bounce happens during a time scale of 
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Figure 3. Slow bounce model (SB). In this model the parameters are: po = —0.586, 
P0 -0.172, px = -0.586, = 2.0, p^x = 0, pxx = 0.5 (with /i = -1, /2 = 1 and 
/a = 1.) Upper panel: background quantities p, (j), Cx and p+p- Lower panel: scale 
factor a and perturbations in the metric ($) for two different wavelengths, fci — 10~^ 
and k2 — 10. 



At ~ 10. The most telling characteristic of the perturbations is that now the mode 
with k2 = 10 experiences more than 20 oscillations during the bounce, while in the MB 
model it only had time to perform about six oscillations. 

The main result of this Section is that cosmological perturbations pass through the 
bounce with their spectrum essentially unchanged. However, our numerical evolution 
cannot address the important question whether there is mixing between the dominant 
and sub-dominant modes, before and after the bounce. This is due to the rapidly 
decaying nature of the sub-dominant solution after the bounce. 

Let us consider the possibility of mode mixing by means of a simplified analytical 

model inspired by the numerically solved fast bounce scenario. Let us take the following 

model for the Hubble parameter: 

H-L -f Hi — H , t , ^ ^, 

H = — + — tanh-, (44) 

2 2 to 

which interpolates smoothly between a de Sitter phase with contraction rate H = H < 
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and an expanding de Sitter phase with expansion rate H = > 0. Notice that 
we have neglected the term ~ k"^, since the numerical analysis have shown that it only 



+ 



+ H+- H_] tanh 
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H+-H^ + {H+ + H^) tanh 
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t ' 
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(45) 



In order to get an analytical solution we assume that = — if_ = h, and we take 
to = 1 for simplicity. (Note that this is strictly equivalent to introducing a new variable 
t/to and rescalling all the constants accordingly.) With these choices, the two linearly 
independent solutions to Eq. IHM are: 



where z 



$1 = 

$2 = 

tanht, P 



(1 

{l-z^f-^^^^ Qliz) 



(46) 
(47) 



v^TTT? > 1, and P;(z) and Q'^Az) are the associated Legendre 



functions of the first and second kind, respectively |32j. 

Since we would like to connect this universe model with a previous contracting 
phase (before pre-cooling) and an ensuing expansion phase (after pre-heating), we should 
consider what happens with the two modes above both at early times {t — * — oo) and at 
late time [t —>■ +oo). This is necessary if we give initial conditions for the perturbations 
and their time derivatives at some initial (early) time, and if we would like to follow 
their evolutions at late times and ask whether that choice of initial values implies a 
mixture of the two modes of Eqs. (l46!) - (H7j) . 

It turns out that the asymptotic limits t — oo {1 + z 0) and t +oo 
[1 — z ^ 0) are very subtle for the associated Legendre functions: in fact, the limit 
1 is singular in the Legendre differential equation, which translates into the fact 
that the Legendre functions are almost degenerate - and, of course, without a complete 
basis of linearly independent functions we cannot accomodate an arbitrary set of initial 
condition. As a result, we need to expand the functions to 3'''^ order in the small 
parameters (1 ± z) ~ exp{—2h\t\) so as to break that degeneracy. The result is that, in 
the limit t — oo {1 + z ^ 0), we get: 

21+/3/2 CSC pTi sin hTT 



X 



(1 + ^) 



r(i-/3) 



r(2 - p) 



l + z 



,2-/3/2 



+ 



$2 ~ 



X 



2i-/3^r(i - p) 



sin hTT T{-h - P)T{1 + h- p)T{l + /?) 
CSC Pit cos hir 



[l + z 



,l+/3/2 



+ 



'l + z 



T{l-P) 

1-/3/2 _ /^(/^ 



i) + 2hr{i-p) 



r(2 - p) 



l + z 



1 2-/3/2 



(4J 



(49) 



K-Bounce 



18 



+ 



1 cot{h + P)n 



Ti-h-P)T{l + h-P)T{l + py ^ ' 



sin hn cos /Svr 

Here we have implicitly assumed that 1 < j3 < 2, so the series in non-integer powers of 
{1 + z) is ordered correctly. 

Comparing the first two terms in Eqs. ( H8l) -( H9l) one can see that they are identical, 
up to a constant. Only the third-order terms in these solutions have different factors. 
Therefore, if we want to assign arbitrary initial conditions to the perturbations at very 
early times, we need to go to third order in the series around {1 + z). Since this small 
factor goes exponentially to zero as exp[A{2 — (3 /2)hti], this means that numerically it is 
very hard to select only one of the modes. Any choice of initial conditions that selected 
one mode at the expense of the other, if made at a very early time ti ^ —^/h, would 
imply a fine-tuning of order exp[4(2 — l3/2)hti]. This means that quite generically, any 
natural choice of initial conditions at very early times will necessarily select a mixture 
of the two modes, and <I>^. 

Consider now the limit t +oo {1 — z ^ 0): 

21+/3/2 

_ (3(p-l) + 2hTil-P) ^^ _ 

4 r(2-/?)^ ' 

^^-^^^f^x{(l-.)-^/^ (51) 
P(P-l) + 2hTil-/3)^^ 2-/5/2 



4 T{2-(3) 



l-z) 



2-^7rcsc/37rr(l + /i + /3) .1+^/2 , 
T{l + h-P)T{l + P) ^ ^' 

It can be noticed also in this limit that the two modes are degenerate up to second order 
in (1 — z). 

The conclusion we can draw from the solutions above at both asymptotic limits is 
that, in order to obtain a pure mode at t ^ 00 one would need to fine-tune the initial 
conditions to order exp[— 4(2 — P/2)h\t\], and inspect the final solution up to the same 
order and precision. It should be evident that a numerical calculation would need an 
astonishing level of accuracy to be able to detect such minute differences. This explains 
why we could not address the question of mode mixing in the numerical analysis. 

With these solutions at hand, one can ask the question: is it actually possible to 
avoid mixing when making a transition between two contracting or expanding de Sitter 
phases? We see on inspection of Eqs. (HHl) to (EH), expanding the functions 1 ± z in 
time, that the general solution for the gravitational potential reads 

^(±) ^ ^(±)e°±* + 5(±)e°±*, (52) 
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where the coefficients A^^\ A'^^\ and a± can be obtained formally from Eqs. fl48p to 
f lSTj) . In the absence of mode mixing, one would expect the transition matrix relating 
i?*^"^)} to \^A^~\B'^~^^ to be diagonal. As the coefficients of the modes are 
different from one side to the other, this is clearly not the case, so one expects mixing, 
at least in this simplified model. 

As a result, if the contracting and expanding phases are of comparable durations, 
an initial condition in the contraction era which is a mixed state of the dominant and 
subdominant modes will become a very different mix of the dominant and subdominant 
modes at the end of the contraction era, in effect transferring power from one component 
to the other. This transfer can result in a large amplification of $. This would mean 
that, provided there is a scale-invariant spectrum in $ before the bounce, however small 
its amplitude, and even if it is present only in the sub-dominant mode, the bounce can 
manage to amplify it to large values. However, this mechanism does not separate the 
spectrum from the amplitude, as, say, the curvaton models [33], because the curvature 
perturbation is, on large, cosmologically relevant scales, conserved through the bounce, 
and thus retains its amplitude as well as its spectral index. 



6. Spectrum of perturbations in a K-bounce 

There are two ways in which we can address the question about cosmological 
perturbations in K-bounce models. First, we could assume that there were no 
perturbations initially, and that a spectrum of cosmological perturbations was generated 
by the bounce itself, through the usual quantum mechanism. Second, we could equally 
well assume that the bounce only distorts a pre-existing spectrum of cosmological 
perturbations. Of course, in general both processes will occur, but in linear theory 
they can be treated separately and the final spectrum will be a combination of the two 
spectra. Let us briefly discuss the first possibility of producing the perturbations at the 
bounce itself. 

While tempting to produce perturbations close to the bounce, one immediately 
faces a major difficulty, namely that it seems rather unlikely that natural, vacuum-like, 
initial conditions could be imposed close to the bounce. Indeed, with the expansion 



( |29|) . one has z = ay —2H / H"^ in ( l20l) . so that switching back to the time variable t, 
we obtain 

where we have set, for simplicity, c\ to unity, since we have seen above (numerically) 
that this quantity is completely regular through the bounce and actually hardly varies 
at all. To leading order in t, Eq. fl2U]) becomes 

V + H^tv + j^— - - J = 0, (54) 

whose general solution is expressible in terms of hypergeometric functions. For all but 
the smallest wavelength modes, these happen to have no oscillatory part to which one 
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could connect the vacuum initial condition fini = e~**^''/v^2fc. The evolution through 
the bounce itself is therefore completely arbitrary for all scales of cosmological interest. 
Hence, in what follows we shall be concerned exclusively with the second possibility, 
namely, of modifying a spectrum originally produced during the contracting phase prior 
to the bounce. 

Recently, |15] considered a K-essence model in which the contracting phase prior 
to the bounce was described by if ~ p/t near the bounce, which leads to rig — 1 2p. 
However, that approximation clearly breaks down at the bounce. We will now assume, 
in the same fashion, that a contracting phase not described by our K-bounce model 
has already taken place, and that a radiation-dominated era starts shortly after the 
bounce. The situation is summarized on Fig. |H We will simply assume that the 
curvature perturbations are transmitted in a non-singular way through these precooling 
and preheating phases - as happens in the usual mechanism of preheating. Therefore, 
here we are just interested in the way an initial spectrum created in a pre-bounce era 
(before precooling) is affected by the bounce, and how that spectrum is transmitted to 
the radiation-dominated era after preheating. Given that all modes of interest are in 
the infrared limit (long wavelengths) before and after the bounce, we only need to ask 
what spectral distortions (if any) are introduced by the bounce. 

The overall conclusion we can draw from the last section is that an initial 
spectrum of long- wavelength modes is entirely unaffected by the bounce: the long- 
wavelength modes all behave in exactly the same way, so their relative amplitudes 
remain unchanged. The curvature perturbation C, is conserved across the bounce for 
large enough wavelengths. What emerges from the numerical analysis discussed in the 
previous section is that for long, i.e. cosmologically relevant, wavelengths, the spectrum 
produced before the bounce, during the contracting phase, is essentially unchanged apart 
from an overall amplification factor for $ which, however, still leaves the curvature 
perturbation constant. Let us consider for instance the case of Fig. [21 very close to 
the bounce, and suppose the contracting phase lasted much longer in the past than 
indicated. Since we are working with arbitrary units in time, we are free to assume that 
the preheating-like phase begins, in this model, around t ~ 7, i.e. very shortly after the 
bounce. Note that we must assume that the expanding era is sufficiently close to the 
bounce in order to ensure a natural solution to the flatness problem. 

Connecting the bounce with the radiation-dominated phase, in a way yet to 
elucidate, very shortly after this bounce took place, we end up with a spectrum of 
perturbations which is undistorted. (recall that our initial conditions were such that 
•I'ini = 1 and $ini = 0.) In other words, the transfer matrix [12], i.e. the matrix that 
relates the initial amplitudes of the growing and decaying modes to their final amplitudes 
at some fiducial instant of time, does not depend on scale. As we have evolved the 
perturbations numerically, it is extremely difficult to extract information about the 
decaying mode. Instead, we focus on the spectrum of the perturbations at some point 
after the bounce has occurred. 

There are two ways in which we can compute the spectrum. First, we can use 
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Figure 4. Typical embedding of the K-bounce into a complete cosmological 
model. To begin with, the Universe is very large, empty, and contracting. At 
that time, the cosmologically relevant modes are all below their potential, i.e. in 
practise their wavelength is smaller than the Hubble scale, and a primordial 
spectrum of perturbation is produced (region marked "PPP" , standing for "Primordial 
Perturbation Production"). Then comes the precooling, when the initial era 
condensates into the effective non-canonical field (j), which then starts to describe the 
cosmological dynamics. A contracting de Sitter stage follows, then the bounce and 
after that the Universe expands in a de Sitter phase. This (presumably short) de Sitter 
era ends through a preheating mechanism, giving way to the radiation-dominated era 



the dominant solution of Eq. ( l26l) . and calculate its amplitude as a function of the 
wavenumber k. By setting all modes to the same initial value, we thus obtain the spectral 
distortions Aj^ caused by the bounce. This is shown in Fig. [5] for the concrete models 
we considered. The main result is that for long-wavelengths (small k) the amplitudes 
are completely flat, meaning that the bounce does not distort the initial spectra. Notice 
also that, for the small wavelengths, the onset of the oscillatory regime influences their 
relative amplitudes, and the change in their spectrum seems to be mo del- dependent. In 
general, slow bounces seem to produce a decay in the spectrum for small wavelengths, 
whereas fast bounces have little to no overall effect over the UV sector of the spectrum 
apart from some oscillations. 

The second method is to evaluate the amplitude of the long-wavelength modes after 
the bounce. This is useful in order to compute the amplification factor A - which is 
only meaningful for those long wavelengths. We have obtained, for the three models we 
studied, ApB =^ 3 x 10^ Amb ^ 3 x 10*' and Age ^ 10^ for the FB, MB and SB models 
respectively - see Figs. 1-3. Note that these values are very much dependent on the 
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Figure 5. Spectral distortions for the FB, MB and SB models. The onset of the 
oscillatory regime is roughly at fc ~ 0.8 for the FB and SB models, whereas for the 
MB model it lies at fc ~ 1.5. The overall normalization in this figure is arbitrary ~ the 
curves only express a change in the relative amplitudes of the modes. 



initial time we put the initial conditions on. This means that in practice, much larger 
amplifications could easily be achieved. However, notice that the physically relevant 
function is the curvature perturbation (. We have found, in agreement with [13], that 
( in fact remains essentially constant - even if $ can be vastly amplified. This means 
that the physical observables are unaffected by the bounce. 

7. Discussions and conclusions 

Bouncing models have been proposed as possible alternatives to infiation. Even though 
such models seem to be able to solve at least some cosmo logical puzzles such as the 
horizon and fiatness problems, many of them still face a basic difficulty of producing an 
almost scale-invariant spectrum of perturbation. The main reason for this failure is that 
there is no generally agreed upon way of making a bounce. In particular, this stems 
from the fact that General Relativity forbids a bounce to take place without spatial 
curvature or violations of the energy conditions. 

Bouncing models have been built based either on a positive spatial curvature [121 [13] 
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or many fluids (TU] - one of tliem liaving negative energy. Tlie results that have 
been obtained up to now are very strongly model-dependent, ranging from no 
effect whatsoever (the modes passing unaltered through the bounce), to a complete 
modification of the spectrum involving fc— mode mixing. It has been argued that 
both the mode mixing and/or the spectral modifications could be due to either the 
spatial curvature and/or the presence of many degrees of freedom, and hence of entropy 
perturbations. Therefore, there has been no agreement about whether bounces would 
lead inevitably to spectral distortions. Hence the need to try and find a simple bouncing 
model with only one degree of freedom and no spatial curvature in the framework of 
GR. 

We have achieved the construction of one degree of freedom simple bouncing models 
by means of a generalized scalar field theory. Our models violate most energy conditions, 
but still we have p > 0. However, because the phantom barrier {p + p < 0) cannot 
be bypassed without some sort of singularity, such models cannot be very reahstic, 
and must be embedded into a more complete theory containing at least the usual 
expanding radiation-dominated era. As we have shown, all our cosmological models 
flow to asymptotic de Sitter solutions with H > 0. This implies that the connection 
with the radiation era must be realized through some sort of preheating mechanism. 
Similarly, the contracting de Sitter solution is a repulsor from which all contracting 
solutions flow. This means that, going backwards in time, all pre-bounce solutions must 
have initiated from a contracting de Sitter stage. Again, in order to relate the bounce 
to a contracting universe dominated by a regular fluid such as dust or radiation, one 
must have a mechanism similar to preheating, but going the other way around, that we 
have called precooling - see Fig. HI 

We have found that, in these simple models, the propagation of perturbations 
is highly non-trivial: although the transition matrix which relates the growing and 
decaying modes before and after the bounce is wavelength-independent (no A;— mode 
mixing in the terminology of the first of Refs. [I2]), it is however non diagonal, so there 
is in general some amount of mixing between the two modes. Analytical calculations 
show that even an exponentially small initial contribution of the sub-dominant mode can 
lead to a high degree of mode mixing in the final spectrum of perturbations. Hence, if 
the two modes have different spectral tilts in the contraction era, the resulting spectrum 
will almost surely consist of a superposition of the two spectra. 

We have also found that this mixing can lead to an amplification of the metric 
perturbations $: any initial suppression of the sub-dominant mode deep in the 
contraction era would be offset by an equal amount of growth prior to the bounce, 
leading to potentially large amplification factors for $. However, this does not impact 
the physically relevant curvature perturbation (, which remains essentially constant 
despite the growth of $, implying that the bounce does not affect physical observables. 
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